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Q"^ I Disks of collisionless particles are important models for certain galaxies and accre- 

■ tion disks in astrophysics. We present here a solution to the stationary axisymmetric 
O , Einstein equations which represents an infinitesimally thin dust disk consisting of two 

streams of particles circulating with constant angular velocity in opposite directions, 

^ijrj' These streams have the same density distribution but their relative density may vary 

• • . continuously. In the limit of only one component of dust, we get the solution for the 

• 5^ \ rigidly rotating dust disk previously given by Neugebauer and Meinel, in the limit 

^ ' of identical densities, the static disk of Morgan and Morgan is obtained. We discuss 

■ the Newtonian and the ultrarelativistic limit, the occurrence of ergospheres, and the 
regularity of the solution. 
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In astrophysics thin disks of colhsionless matter, so called dust, are discussed as models 
for certain galaxies (see e.g. [|l|) or for accretion disks. A fully relativistic treatment of 
these models is necessary if a black hole is present since black holes are genuinely relativistic 
objects. But also the exact treatment of dust disks without central object would provide 
deep insight both in the mathematical structure of the field equations and in the physics 
of rapidly rotating relativistic bodies since dust disks can be viewed as a limiting case for 
extended matter sources. Corresponding exact solutions hold globally - in the vacuum and 
in the matter region - and can thus provide physically realistic testbeds for numerical codes. 
Since Newtonian dust disks are known to be unstable and since there are hints by numerical 
work (see e.g. 0) that the same holds in the relativistic case, such solutions could be taken 
as exact initial data for numerical collapse calculations. Whereas the Newtonian theory of 
such disks is well established (see [0] and references given therein), the same holds in the 
relativistic case only for static disks which can be interpreted as consisting of two counter- 
rotating streams of matter with vanishing total angular momentum. The first disk of this 
type was considered by Morgan and Morgan 0. Infinitely extended dust disks with finite 
mass were studied by Bicak, Lynden-Bell and Katz in the static case and by Bicak and 
Ledvinka [^] in the stationary case. The first to construct the exact solution for a finite 
stationary dust disk were Neugebauer and Meinel who gave the solution for the rigidly 
rotating dust disk which was first treated numerically by Bardeen and Wagoner 0. They 
solved the corresponding boundary value problem for the Einstein equations with the help 
of a corotating coordinate system. 

In this letter we present a class of new disk solutions to the Ernst equation where such 
a coordinate system cannot be used. The disks of finite radius po consist of two counter- 
rotating components of dust with respective density a'^{p). The angular velocity of both 
streams of particles is of the same constant absolute value VL but of different sign. The 
relative density 7 = (o"+ — a^)/ (cr+ + o" ) is a constant with respect to p and C which varies 
between one, the rigidly rotating dust disk |]^, and zero, the static Morgan and Morgan 
disk 0. Interestingly, observations indicate that the galaxy NGC 4550 is built from 
two counter-streaming stellar components. We are able to give the explicit solution for the 
above configuration in dependence of the three parameters po, ^ and 7 which parametrize a 
Riemann surface of genus 2. We discuss physically interesting features like the static limit, 
the Newtonian regime and the ultrarelativistic limit. We study the occurrence of ergospheres 
and the absence of singularities in the allowed range of the physical parameters. 
Newtonian dust disks 

It is instructive to consider first the Newtonian case where the gravitational field is given 
by a scalar potential U which is a solution to the Laplace equation Af/ = 0. The potential 
U has to be everywhere regular except at the disk where the balancing of the centrifugal 
and the gravitational force leads to boundary values for U . We use dimensionless (po = 1) 
cylindrical coordinates (p, C, <P) and put the disk in the C, = 0-plane which leads at the disk 
to 

u,ip,o) = n'p. (1) 

For constant fl, this can be easily integrated to give U (p) = Uo + |^^P^ where the constant 
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Uq is related in the relativistic case to the central redshift. Notice that there are no effects 
due to counter-rotation in the Newtonian case since fl enters (0) quadratically. The solution 
to this problem can be constructed e.g. with Riemann-Hilbert techniques (see ^) which lead 
in the case of the Laplace equation to a potential of the form 

4iri Jt ^(r - C)2 + 

which is a function on the Riemann surface of genus zero given by hI{t) = (t — C)^ + p^; here 
In G{t) is an analytic function and F is the covering of the imaginary axis in the upper sheet 
between — i and i. In the Newtonian case, the situation is equatorially symmetric which 
leads at the disk to 

1 np InGMdr , , 

/vri Jo vr^ + 

This is an Abelian integral equation for In G since the left hand side is given by the integral 
of (|l|). It has the solution InG = 4^2^ (r^ + 1) where a singular ring at the rim of the disk 
is excluded by the condition ^(±1) = 1. The condition implies that U is continuous at the 
rim, and it fixes the constant f/o = This shows that the solution indeed only depends 

on the two parameters po and fi. 
The relativistic case 

The metric describing the exterior (i.e. the vacuum region) of an axisymmetric, stationary 
rotating body can be written in the Weyl-Lewis-Papapetrou form (see |^) 

ds^ = -e^^idt + ad0)2 + e-2^ {e^\dp^ + dC) + pW) , (4) 

where dt and are two commuting asymptotically timelike respectively spacelike Killing 
vectors. In this case, the vacuum field equations are equivalent to the Ernst equation for the 
potential /, 

1 2 

where / = + ih and the real function h is related to the metric functions via hz = 
— (i/ p)e^^ ttz- Here the complex variable z stands for z = p + iC,. With a solution of the Ernst 
equation, the metric function U follows directly from the definition of the Ernst potential 
whereas a and k can be obtained from / via quadratures. 

The complete integrability of this equation allows in principle to give explicit solutions 
for boundary value problems. For two-dimensionally extended matter sources, e.g. infinitesi- 
mally thin disks, one gets global solutions since one encounters ordinary differential equations 
in the matter region which provide boundary data for the vacuum equations. Riemann- 
Hilbert problems (see e.g. 0]) generate solutions to the Ernst equation with free functions 
which have to be determined by the boundary data. However this is only possible uniquely 
in the case of Cauchy data, which would lead in general to singular solutions of the elliptic 
field equations. In p[ this problem was circumvented by the use of a corotating coordinate 
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Figure 1: The homology basis for E2. 



system, a technique that is not apphcable in the case of differential rotation or for several 
matter components. 

In we have shown that the matrix Riemann-Hilbert problem is gauge equivalent to 
a scalar problem on a Riemann surface which can be explicitly solved in the case of rational 
boundary data, i.e. belongs to Korotkin's ||T^ finite gap solutions. The explicit form of the 



solution allows in principle to solve boundary value problems directly. The branch points of 
the Riemann surface are used as a discrete degree of freedom. Since the solutions already have 



the required regularity properties (see []12| , pT^), there can be no problems with singularities 
as in the direct approach in the matrix case. The purely algebro-geometric treatment has 
the further advantage that it is not limited by the possibility to introduce special coordinate 
systems. 

The simplest surface leading to equatorially symmetric Ernst potentials is of genus 2 and 
of the form iy^{K) = {K - Pq){K - Po){K^ - Ef){K^ - Ef) where Pq = -iz is a moving 
branch point related to the physical coordinates whereas the ^-independent branch points 
are given by Ei = — (ai + i/3i) with positive ai and Pi. We introduce the standard quantities 
associated with a Riemann surface, with the cut system of figure ^ namely the 2 differentials 
of the first kind doOi normalized by duj = 27ri(5jj, the Abel map uJi{P) = Jp^ dui which is 
defined uniquely up to periods, the Riemann matrix 11 with the elements Hij = dcUj, and 

the theta function 9(2;) = Sivgz^ ^xp || (HA^, A^) + {z, A^)|. The normalized (all a-periods 
zero) differentials of the third kind with poles in p and q and residues +1 respectively —1 
are denoted by dujpq. Then one can show that 

/(p.O- ^i"!°°!'t"! exp(-L /lnGMd.^.„-M|, (6) 

B(t<;(oo ) + u) y zni J J 

where the two-dimensional vector u has the components = 2^ /r 1^ Gduji and F as in the 
Newtonian case, is a solution to the Ernst equation which is everywhere regular except at 
the disk if 6(u;(cx)~) +u) 7^ 0. We can also give the metric function a (see and |]l3l) and 
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k (see [n| and references given therein) explicitly. Notice that solutions of the form are 



characterized by one real valued function G and two real numbers, e.g. a and P defined by 
Ef =: a + ip where /? has to be positive. 

Dust disks within general relativity are best described with the help of Israel's covariant 



formalism where the regions C > ^ind C < are matched at the hypersurface C = 0- The 
surface stress-energy tensor of the disk, for counter-rotating dust Sab = <^^v^v^ + ^'""'^a ^(T 
where vf is the velocity of the respective component of the the two-dimensional dust, is 
related to the jump of the extrinsic curvature at the disk. If we put Z = {a + 7/^)6^'^, 
A = 2Q'^e~^^° and 6 = [1 — this leads via the Einstein equations to the boundary 

conditions 

^C" 2^ ^^-^ 2^ 

The second equation may be integrated to give Z"^ — + 5e^^ = ^e^^. Notice that A can 
be viewed as a parameter that indicates how relativistic the situation is: the larger it is, the 
larger is the angular velocity and the central redshift. 

We can now state the following theorem which is the main result of this letter. 
Theorem: 

The solution of the boundary value problem for the Ernst equation is given by an Ernst 
potential of the form (P) with 

and 



G{r) = . (9) 



We briefly sketch the proof which uses in principle only fundamental theorems of algebraic 



geometry which can be found in the standard literature (see e.g. ||T^, [T^)- The main idea 
is to establish the relations between the real and the imaginary part of the Ernst potential 
enforced by the given Riemann surface independently of This selects the class of bound- 
ary value problems that can be solved on a certain Riemann surface. We note that the steps 
below can all be performed uniquely what implies that the form of the solution is unique at 
least among genus 2 solutions. 
Proof: 

1. The Jacobi inversion theorem ensures that the equations uj{X) — oj{D) = u where 
X = Xi + X2 and D = Ei + {~Ei) are divisors can always be solved for the divisor 
X. The divisor X is called non-special if the solution is unique (how to overcome problems 
with special divisors is discussed in [|T^). If we introduce this divisor in (^), we end up with 

^The simplest example for this concept is provided by the function e"" which is the analogue of theta 
functions on the Riemann surface of genus zero. In this case we have of course the simple relation for the 
real and imaginary part that cos^ w + sin^ w = 1 for all w. 



4 



the algebraic formulation of the hyperelliptic solutions in |T8| and |jT9[. 

2. The reality of the Ui leads to the condition uj{X) +uj{D) = uj{X) +uj{D). Abel's theorem 
implies that this condition is equivalent to the existence of a rational function on S2 with 
zeros in X + D and poles in X + D. We can thus express X via he'"^^ and a second real 
quantity as the solution of a system of algebraic equations. 

3. If we introduce the divisor X in the expression for the metric function a ([11, 
use so-called root functions (see |T^) to express X via a and 6e~^^ alone. 

4. Differentiating with respect to the physical coordinates and using the equatorial symme- 
try in the equatorial plane, we arrive at a set of equations which contain only e^^, a, h and 
their derivatives. Using the boundary conditions (^ together with the definition of 6, we 
can show that these equations are identically satisfied. 

5. The constant (3 is used as an integration constant in the integration of the second equation 
in(0). 

6. The function G is determined as in the Newtonian case as the solution of an Abelian 
integral equation, e.g. the equation following from uJi{X) — ui{D) = ui. The regularity 
condition G{±i) = 1 then fixes a. 



Discussion: 

The above disk solutions have several interesting limiting cases. The simplest is of course 
6 = 0, the rigidly rotating dust disk §], here in the form |jl2|. Mathematically more intriguing 



is the limit of the Morgan and Morgan disk. It was argued in |]T2[ that the hyperelliptic 
solutions of the form (|^) are all non-static. This holds of course only in the case of regular 
Riemann surfaces. A real Ernst potential, i.e. a solution which belongs to the static Weyl 
class, is only possible on a degenerated Riemann surface where the branch points coincide. 
This is a well-known limiting case of algebro-geometric solutions of integrable equations in 
which the soliton solutions are obtained. Since the solutions via Backlund transformations 
correspond to the 'solitonic' solutions of the Ernst equation, Korotkin [|TT| was able to obtain 
the Kerr solution in such a limit. 

The problem of the parametrisation of the solution via A and 6 which is enforced by the 
solution process is that the dependence on the physical parameters fl and 7 is concealed by 
a factor e'^^° which is itself a hyperelliptic function of 6 and A. In the static limit, we get 



however the simple relation 6 = 2(1 + y 1 + 1/A^) where (3 = 0. This means that the branch 
points coincide pairwise on the real axis since a > in this case. The Ernst potential thus 
simplifies to 

l/(p.C) = -J-/ /■'g'^'''^ (10) 

4m Jt ^(r - Cr- + p' 

where G = 1 — 4(r^ + l)/(5, i.e. the expected potential of the Morgan and Morgan disk in the 



form ||2^. This disk thus provides a nice example where a stationary solution of the Einstein 
equations for an extended body is linked continuously to a static solution via a parameter. 

The Newtonian limit of the disks is reached for A ^ 0. In this case, we have G ~ 1 
which implies that the solution approaches Minkowski spacetime as expected. Since In G ~ 
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Figure 2: Ergospheres in the disk in dependence of A and 6. 



2A(r^ + 1) and U is given by (|^), we get the above Newtonian hmit. As expected, counter- 
rotation i.e. 6 does not play a role in this case. The other extreme limit is the ultrarelativistic 
limit. The corresponding limit of the Morgan and Morgan disk is of course static. It is 
reached for 5 — 4 or A ^ oo and describes a disk where the matter streams rotate with 



the velocity of light (see e.g. |23). The density of the matter diverges at the origin where 



the central redshift diverges, too, but the mass of the disk remains finite. For the non-static 



disks, the exterior of the disk can be interpreted as the extreme Kerr solution as in @, p^ . 
The disk can thus be viewed as being hidden behind the horizon of the Kerr metric. The 
exact form of these extreme configurations will be the subject of further investigation. 

In JT^ , |T^ we have given the necessary conditions for the occurrence of ergospheres. 
Numerically we find that the ergospheres always form first in the disk. If more matter is 
counter-rotating, ergospheres are less likely to appear since they are due to gravitomagnetic 
effects which vanish in the static limit of two counter-rotating streams with identical density. 
The necessary condition for an ergosphere to reach the rim of the disk is that A takes the 
value 2/(1 — S) which is only possible for 6 < 1. For larger values of 6, ergospheres thus 
cannot meet the equatorial plane outside the disk. In figure 2 we have plotted for given 6 
the values A, p of the common points of the ergosphere and the disk. For small values of A, 
there will be no ergospheres since we are close to the Newtonian regime. For increasing A, 
there will form a ring-like ergosphere at some radius p. If A is increased further, the then 
roughly toroidal ergosphere hits the disk at two values of p. The larger value of p reaches 
the rim of the disk at a finite value of A for 5 < or is strictly smaller than 1 for 6 > 1. The 
inner radius of the ergosphere will hit the axis in the ultrarelativistic limit for A = Ac. 

It makes little sense to increase the parameter A beyond this critical value Ac where the 
solution still formally exists. These regions are most probably non-physical since the solu- 
tions may have negative ADM-mass (see [|12|) and have singularities. The allowed parameter 
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range is thus < A < Ac and < S < 2(1 + y 1 + l/A^). It can be shown that the regularity 
condition Q{uj{oo^) + n) 7^ is always fulfilled in the allowed parameter range. 

Thus we have shown that it is possible by purely algebro-geometric methods to construct 
a solution to the Ernst equation which describes a non-static counter-rotating dust disk. 
The fact that this is possible without the use of special coordinate systems gives rise to the 
hope that astrophysically interesting objects like dust disks with differential rotation or with 
a central black hole or disks with surface tensions can be constructed in selected cases. The 
open question is whether the resulting algebraic equations can still be handled analytically. 
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